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Abstract 

Boundary integral equation is derived for the problem of scattering 
of electromagnetic waves by 3D homogeneous body of arbitrary shape. 
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1 Introduction 

Let D C be a bounded domain with a smooth connected boundary S. 

The exterior domain D' = \ D is filled with a homogeneous material 

with parameters (eo,cr = 0,/io), D is filled with a material with parameters 

C g _l_ i£ in 

(e,(T, /xo), e, eO) cTj A'o are constants. Let e' = < ^ ■ i^'i Here u 

[ eo, m V . 

is the frequency, no is magnetic parameter, e, eo are dielectric parameters, 
(T > is conductivity. 

The governing equations in are 

V X E = iunoH, V X H = -iuje'E. (1) 

The boundary conditions on S are continuity of the tangential components 
of E and normal components of e'E across S: 

[N,E+] = [N,E-], (2) 
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Euclidean basis, 5ij 



^ k^inD\ fc2-' .2..„/ (6) 



and 

N-{e'E+) = N-{eoE-). (3) 

Here N is the unit normal to S, pointing into D', [N, E] (E ■ N) is the 
cross (dot) product of vectors, and E^{E^) is the limiting value of E as 
x^s£S,xe D{D'). 

We assume that the incident field {Eq,Hq) satisfies equations ([T]) in 
with e' replaced by eo- For example one may take the incident field to be 
a plane wave: Eq = eie*'^^^, where Cj • ej = 6ij, {ej}j^i is the standard 

1, if ^ = J, 
0, if i j. 
If E is found then 

H = {iojfior^^ X E. (4) 

From ([T|) one gets 

V X V X E - JC'^E = in R^, (5) 

where 

in D, K^ = u'^e'no, 

The scattering problem consists of finding the solution of equations (l5|) , ([2]) , 
([3]), such that 

E = Eo + V, (7) 
Vr - ikV = o (^^^ , r:=|a;|^cx3. (8) 
Assumption ([8]) means that V satisfies the radiation condition. 

2 Derivation of the boundary integral equations 

Let us look for the solution to (P, §1, of the form: 

( V X JsGix,t)J{t)dt, xGD 

\ V X fgg{x,t)j{t)dt + Eo{x), X G D', ^""^ 

where 

„ik\x—y\ „iK\x—y\ 

g{x,y) = — G{x,y) = — (10) 

The currents j and J are vector fields tangential to S. 

Thus, there are four scalar unknowns: the two unknown vector fields j 
and J, tangential to S. 
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For any j and J vector E solves equation (l5|) in D and in D' and satisfies 
conditions ^ and ([8]) because g satisfies the radiation condition ([8]). Thus, 
([9|) is the solution to the scattering problem if j and J can be chosen so that 
the boundary conditions ([2]), ([3]) are satisfied. 

Condition ([5]) is a vector equation, which is equivalent to three scalar 
equations, and equation ([3]) is a scalar equation. Therefore equations Q 
and (l3|) together are equivalent to four scalar equations for four unknown 
scalar functions, the coordinates of the tangential to S vector fields j and 
J. 

Equation ([2]) can be written as a Fredholm-type integral equation. We 
have 



N^,[V^G{x,t),J]]\^^s, xeD = -^-4^^+ [ VsG{s,t)N,- J{t)dt, 

^ Js 



(11) 

where the known formula for the limiting value of the normal derivative of 
the single-layer potential was used (see, e.g., f2]). 

Since J{s) ■ Ng = (because J(s) is a tangential vector field) and J{s) is 
assumed Lipschitz, the last integral in (llip converges absolutely. Since we 
assume the surface sufficiently smooth, e.g., S G C^'" a > 0, and the incident 
field is smooth, the currents j and J are as smooth as the data, in particular, 
they are Lipschitz. The class of surfaces, satisfying the condition S € C^'", 
consists of surfaces whose graph in local coordinates is differentiable and its 
derivative satisfies the Holder condition with the exponent a £ (0, 1]. 

The operator in (jlip is defined as 

Thus, equation ([2]) can be written as: 
_^4+J + J^ / VMs,t)N,-J{t)dt = -^^^ 

J S 

+ f Vsgis,t)Ns-jit)dt + [N,Eo], 
Js 

(13) 

where 
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Equation (jl3p is of Fredholm type: the integral operators in ()13p are compact 
in C{S) and in L?'{S). Equation ([3]) yields: 

Ns- [ [V^G(x,t), J(t)]U^,, ^eodt = N^- [ [Vxg{x,t),j{t)]\^^s, xeD'dt+Ns-Eo. 
Js Js 

(15) 

Equation (fTSj) is singular. 

Claim: The integrals in (|15|) exisi as Cauchy principal values. 
Let us verify this claim. One has 



(16) 



Ns ■ [VMx,t),Jit)]^^D, xGD = Ns ■ [^-^{iK Vts,J{s)] 

47rrst rst 

+ O ( — ) , := — , Tst = \rtsl 
\rstj Tst 

because \ J{s) — J{t)\ < c\s — t\. 
The singular term in (jl6p is 

-i^iV^.[rO„J(s)] = ^sin0, (17) 

where 6 = 0{s,t) is the angle between the x-axis and the vector fgt- We 
choose the x- axis in the plane tangential to S at the point s so that it is 
directed along the vector J{s). 
Since 

sin edd) = 0, t = e''^, 0<6<2tt, 







the Claim follows from Theorem 1.1 on p. 221 in [T]. This theorem says 

\x-y\^' 



that a singular integral /Jf'^,l, u{y)dy exists as a Cauchy principal value 



if 

y -X 



f{x,e)ds = 0, 



ism-i \y-x\ 

Numerical methods for solving Fredholm equations and singular integral 
equations are well developed ([!])• They are not discussed here. 
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